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This  report  summarizes  the  result  of  the  research  carried  out 
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study  of  the  propagation  characteristics  of  wave  passing  through  an 
ensemble  of  nonspherical  particles  for  which  the  multiple  scattering 
effects  are  important.  Using  the  vector  transport  theory  approach 
which  includes  the  depolarization  effects  due  to  irregularly  shaped 
particles,  the  incoherent  specific  intensities  caused  by  the  multiple 
scattering  effects  were  computed.  Also  described  in  this  report  are 
several  future  research  areas. 
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I.  Introduction 

This  is  a  final  report  on  a  study  sponsored  by  the  Army  Research  Office 
(DAAG  29-82-C-0008)  from  April  12,  1982  through  April  11,  1984.  The  principal 
objective  of  this  research  was  to  learn  the  propagation  behavior  of  wave  passing 
through  an  ensemble  of  nonspherical  particles  for  which  the  multiple  scattering 
effects  are  important.  Using  the  vector  transport  theory  approach  which  includes 
the  depolarization  effects  due  to  irregularly  shaped  particles,  we  were  able  to 
obtain  the  incoherent  specific  intensities  caused  by  the  multiple  scattering 
effects.  This  final  report  gives  a  summary  of  our  accomplishments  during  this  phase 
of  the  research  program. 


II.  Summary  of  Accomplishments 

Results  of  our  investigation  are  summarized  in  the  following: 

(1)  Matrix  representations  of  the  vector  radiative  transfer  theory  for  randomly 

2 

distributed  nonspherical  particles. 

In  recent  years  considerable  attention  has  been  paid  to  the  problem  of  communica¬ 
tion  through  various  particles  in  the  atmosphere  such  as  fog,  ice  particles,  snow, 
etc. **  These  particles  are  often  nonspherical  causing  polarization  dependent  attenua¬ 
tion,  depolarization,  and  cross  polarization.  These  polarization  effects  are  also 

7_8 

evident  in  lidar  and  radar  dectection  and  terrain  scattering. 

In  order  to  include  complete  polarization  effects  in  the  radiative  transfer  theory, 

the  equation  of  transfer  must  be  expressed  in  matrix  form  using  the  Stokes'  vector. 

5-9 

Even  though  the  vector  radiative  transfer  equation  has  been  studied  previously,  most 
of  these  studies  are  directed  to  spherical  particles  resuling  in  considerable  simpli¬ 
fication  of  the  equation.  For  example,  the  extinction  matrix  becomes  a  scalar,  and 
the  Mueller  matrix  becomes  a  function  of  <p-4>’.  For  nonspherical  particles,  the 
coherent  field  becomes  depolarized.^  We  have  derived  a  compact  matrix  representation  of 
the  propagation  of  the  coherent  field  in  terms  of  the  Stokes'  vector  as  well  as  a  compact 


representation  of  the  extinction  and  Mueller  matrices  using  the  direct 
product.  These  formal  matrix  representations  may  facilitate  calculation  of  the 
propagation  and  scattering  characteristics  of  the  Stokes*  vector  for  nonspherlcal 

particles.  Results  of  this  investigation  are  given  in  a  paper  which  has  been 

2 

submitted  for  publication  in  the  Journal  of  the  Optical  Society  of  America. 

3 

(2)  First  order  multiple  scattering  theory  for  nonspherlcal  particles. 

It  is  known  that  the  atmosphere  contaminated  by  smoke  or  hydrometeors  (ice 

or  rain)  has  profound  effects  on  the  propagation  of  high  frequency  electromagnetic 

energy,  which  in  turn  may  degrade  the  performance  of  many  modern  communication 

1 

or  detection  systems.  Smoke  or  hydrometeors  may  take  on  different  shapes  other 

than  spheres  and  the  close  proxlminity  of  these  particles  may  enhance  multiple 

7,11,12 

scattering  effects.  It  is  therefore  important  to  assess  the  relative  im¬ 

portance  of  nonspherlcal  versus  spherical  particle  effects  in  the  multiple  scat¬ 
tering  of  electromagnetic  waves  by  an  ensemble  of  particles. 

Various  studies  of  the  multiple  scattering  effects  of  hydrometeors  on  micro- 

wave  propagation  have  been  conducted*2-*^  While  earlier  studies  dealt  with 

14 

scalar  fields  and  first-order  solutions,  more  recent  ones  have  considered 
circularly  or  linearly  polarized  incident  waves  and  exact  solutions  to  the 
aquation  of  transfer.  In  all  of  the  previous  studies  only  spherical  particles 

were  considered.  The  emphasis  in  this  investigation  is  on  the  effects  of  non- 
apherlcal  particles  as  compared  *wlth  spherical  particles  on  the  coherent  as  well 
as  incoherent  intensities  produced  by  multiple  scattering. 

Chur  formulation  for  the  incoherent  intensities  is  based  on  the  equation  of 
transfer  for  the  Stokes'  parameters,  while  the  formulation  for  the  coherent  In¬ 
tensities  is  based  on  the  classical  approach  of  van  de  Bulst.  A  linearly  polari¬ 
sed  wave  Is  assumed  to  be  obliquely  incident  upon  a  plane-parallel  medium  con- 


2 


talning  the  ensemble  of  particles.  In  principle,  the  equation  of  transfer  can 
be  solved  exactly  using  the  eigenvalue-eigenvector  technique  with  the  scattering 
characteristics  of  the  particles  being  calculated  using  the  extended  boundary 
condition  method^  and  a  given  drop  sice  and  orientation  distribution,  numerical 
computation  based  on  this  technique  is  being  carried  out  and  results  vlll  be  pre¬ 
sented.  The  main  purpose  of  the  present  Investigation  is  to  present  the  first- 
order  solution  to  this  problem  dealing  with  non-spherical  particles.  It  is  be¬ 
lieved  that  providing  the  first  order  expressions  for  the  incoherent  intensities 
represents  the  necessary  initial  step  towards  the  complete  solution  of  the  multi¬ 
ple  scattering  problem. 

As  specific  examples,  numerical  results  based  on  the  first-order  solution 
are  obtained  for  the  co-polarized  and  cross-polarized  intensities  as  functions 
of  observation  angle  and  optical  depth  for  two  types  of  particles:  those  with 
high  loss  tangent  (l.e.,  high  absorbent)  and  those  with  low  loss  tangent  (l.e., 
low  absorbent).  As  expected,  we  found  that  multiple  scattering  effects  are 
more  pronounced  for  ensemble  of  low  absorbent  particles  than  for  ensemble  of 
high  absorbent  particles. 

Kesults  of  this  investigation  are  given  in  a  paper  which  has  been  submitted 
to  Applied  Optics  for  publication? 

(3)  Vector  radiative  transfer  equation  for  nonspherlcal  particles. 

It  is  important  to  distinguish  the  vector  radiative  transfer  equation 
which  is  capable  of  including  the  polarization  effects,  from  the  scalar  rad¬ 
iative  transfer  equation.  In  the  following  we  shall  delineate  the  formula¬ 
tion  of  the  vector  transfer  equation  for  an  ensemble  of  nonspherlcal  particles. 

Let  us  now  consider  a  specific  geometry  in  order  to  obtain  the  multiple  scat- 
tering  effects.  The  geometry  consists  of  a  layer  of  uniformly  distributed 
nonspherlcal  particles  of  thickness  d.  (See  Fig.  1)  A  linearly  polarized 


coherent  plane  wave  is  incident  upon  this  layer  of  particles.  Ue  are  inter¬ 
ested  in  the  transmitted  co-polarlzed  coherent  intensities  as  well  as  the 
transmitted  co-polarlzed  and  cross-polarized  incoherent  Intensities.  Con¬ 
sider  the  modified  Stokes'  parameters  (1^,  ,  U,  V)^: 

*1  "  'hh'-h  *  <E2E2%  ^  (1) 

U  -  2Re<E1E2*>,  V  -  2  > 


where  Ej  and  E2  are  the  electric  field  components  in  the  e  ♦directions, 
respectively  (Fig.  1),  Re  and  Im  denote  the  real  part  and  the  imaginary  part 
and  the  asterisk  means  the  complex  conjugate.  The  equation  of  transfer  for 
arbitrarily  shaped  particles  is  given  by 

IT)U)  ♦  J IsHi’Hu,'  ♦  1^1  <2> 

where 

U1  - 


which  is  a  4x1  incoherent  specific  intensity  matrix,  and  where 
[I  ]  “  4x1  incident  specific  intensity  matrix 
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and  (8,$)  correspond  to  the  incident  and  ecettered  veve  directions, 
respectively  end  £  is  the  angle  between  the  Incident  and  scattered  waves. 
The  scattering  functions  A|j,,  A|r«  A*^,  and  A^  from  a  single  arbitrarily 
shaped  scatterer  can  be  calculated  according  to  the  extended  boundary 
condition  approach  (the  T-matrix  approach).15 

The  Incident  specific  intensity  Skatrix  [1^]  Is  the  source  function 
generated  by  the  scattering  of  the  reduced  Incident  specific  intensity 
(1^1  propagating  in  the  direction  of  y^  *  cos  6^  and  ■  0,  and  is  given 


ky 
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where  £ 1  . 1  is  given  by  the  incident  Stokes  vector  [I  ] 
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and  (Q]  is  the  4x4  matrix .  Ve  can  also  express  [2rl1  using  the  field 
representation . 


If  the  incident  wave  with  the  total  intensity  of  unity  is  polarised  parallel 
to  the  x-x  plane  (vertical  polarisation),  (IqI  is  given  by 
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Similarly,  if  the  incident  wave  is  pola.lxed  perpendicular  to  the  x-x 
plane  (horizontal  polarization) ,  we  have 


Ho1 


(17) 


The  boundary  condition  for  the  incoherent  specific  intensity  [i]  is 
ll)  ■  0  cos  e  >  0(0  i  V  S  1)  at  t  ■  0 

(l)  •  0  cos  6  *  0(0  iVM)  at  *  -  d 
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(18) 


The  equation  of  transfer  (Eq.(2))  and  the  boundary  condition  (18) 

constitute  the  complete  mathematical  formulation  of  the  problem.  This  is 

an  Integra-differential  equation  with  five  variables  (zt  y,  y',  and  4>')* 

(4)  Multiple  scattering  calculations  for  nonspherical  particles  based  on  the 

4 

Vector  radiative  transfer  theory. 

In  this  investigation  we  shall  consider  the  multiple  scattering  of  a 
linearly  polarized  plane  wave  Incident  upon  a  plene-parallel  slab  consisting 
of  uniformly  distributed  identical  particles.  It  is  assumed  that  every  part¬ 
icle  is  a  spheroid  whose  symmetry  axis  is  normal  to  the  slab.  Ve  aim  at  com¬ 
puting  the  incoherent  field  that  is  generated  within  the  slab,  in  terms  of 

its  Stokes'  vector.  Our  starting-point  is  the  integro-dif f erential  equation 

5.9 

of  radiative  transfer.  ’  This  formulation  provides  the  incoherent  Stokes' 
vector  for  a  given  normal  or  oblique  illumination  as  a  function  of  the  scatter¬ 
ing  amplitudes  of  a  single  particle.  Both  attenuation  and  multiple  scattering 
are  taken  into  account.  The  analysis  is  based  upon  the  Fourier  expansion  of 
this  equation  of  transfer.  Because  of  the  specific  symmetry  and  orientation  of 
the  particles,  the  Fourier  components  of  the  Stokes'  vector  are  independent  of 
one  another.  Such  a  decoupling  enables  us  to  determine  a  given  component  by 
solving  the  separate  equation  of  transfer  that  it  will  satisfy.  After  an 
estimate  of  the  Integral  terms  by  Gauss'  quadrature,  the  component's  values  on 
a  net  of  discrete  backward  and  forward  directions  are  solutions  of  a  system 
of  first-order  differential  equations  with  constant  coefficients  and  are  com¬ 
puted  by  an  eigenvalues-eigenvector  technique.  A  partial  sum  of  the  Fourier 
series  composed  of  such  components  represents  the  Incoherent  Stokes'  vector; 
in  the  case  of  normal  Incidence,  only  tow  such  components  are  different  from 
zero  (order  0  and  2). 

The  behavior  of  the  Incoherent  field  is  thereby  investigated  for  low-loss 
particles  (ice)  and  high-loss  particles  (smoke)  illuminated  by  x-polarized  waves 


under  typical  incidences.  Patterns  of  the  x-polarized  (co-polarized)  and 
y-polarized  (cross-polarized)  incoherent  Intensities  versus  the  direction 
of  observation  are  computed.  These  results  enable  us  to  examine  the  con¬ 
vergence  of  the  process  and  the  Influence  of  the  particles'  geometry  and 
distribution.  The  consequences  of  the  low-  or  high-losses  are  emphasized. 
Comparisions  with  results  obtained  by  a  first-order  scattering  theory 
are  also  carried  out.  Initial  indication  is  that  excellent  qualitative 
as  well  as  quantitative  agreement  was  found  for  low  particle  density  -  small 
optical  depth  medium  and  good  qualitative  (but  not  so  good  quantitative) 
agreement  was  also  apparent  for  higher  particle  density  -  larger  optical 
depth  medium. 

Results  of  this  investivation  will  be  given  in  a  paper  which  will  be  sub- 

4 

mitted  to  Radio  Science  for  publication. 

Let  us  now  discuss  further  the  two  ways  of  solving  the  vector  radiative 
transfer  equation  as  discussed  in  sections  (2)-(4)  above.  The  first  way,  the 
exact  way,  starts  with  the  expansion  of  the  vector  radiative  transfer  equa¬ 
tion  in  Fourier  series  in  The  integral  with  respect  to  y*  is  converted  to 
a  series  representation  by  the  Gauss'  quadrature  formula.  The  resultant  set 
of  coupled  linear  first-order  equations  are  then  solved  by  the  eigenvalue- 
eigenfunction  technique  with  the  given  boundary  conditions.  The  second  way, 
the  approximate  way,  uses  a  perturbation  technique  on  the  transfer  equation. 

The  result,  called  the  first-order  solution,  is  a  set  of  analytic  expressions 
for  the  incoherent  intensities.  We  have  progressed  in  both  of  these  directions 
The  attractiveness  of  the  first-order  approach  is  in  its  simplicity  and  clarity 
and  the  case  of  obtaining  numerical  results  while  the  exact  approach,  al¬ 
though  very  complicated,  provides  accurate  results  which  can  be  used  as  ref- 


erence  or  check  points  for  the  first-order  or  other  approximate  approaches. 

The  exact  results  may  also  provide  regions  of  validity  for  the  approximate 
first  order  expressions.  It  is  noted  that  because  of  complexity  of  the 
exact  approach  and  the  many  resultant  numerical  difficulties  associated  with 
large,  ill-conditioned  matrices,  extreme  care  must  be  taken  in  performing  the 
calculations.  We  have  successfully  treated  the  problem  of  a  slab  of  uniformly 
distributed  Identical  oblate  or  prolate  spheroidal  particles  with  their  sym¬ 
metry  axis  aligned  along  the  normal  direction  of  the  slab.  Different  sizes, 
shapes,  and  densities  of  the  particles  as  well  as  different  incident  angles 
for  the  plane  wave  were  considered.  This  is  the  first  time  that  such  a  cal¬ 
culation  taking  into  consideration  the  complete  polarization  effects  has  been 
carried  out  for  an  ensemble  of  the  nonspherical  particles. 

It  can  be  seen  from  the  above  summary  of  accomplishments  that  starting 
with  the  vector  radiative  transfer  theory  including  the  depolarization  effects 
we  have  been  able  to 

•  obtain  the  first-order  expression  for  the  coherent  and  incoherent 
intensities  for  waves  propagating  through  an  ensemble  of  spherical 
or  nonspherical  particles, 

•  expand  the  vector  radiative  transfer  equation  for  nonspherical  part¬ 
icles  in  Fourier  series.  (In  general,  all  Fourier  components  are 
coupled  to  each  other.  We  have  Investigated  the  special  case  for 
which  the  Fourier  components  are  decoupled), 

•  obtain  specific  numerical  solution  to  the  vector  radiative  transfer 
equation  for  certain  orientation  and  symmetry  of  the  nonspherical 
particles, 


•  discuss  the  limitation  of  the  radiative  transfer  theory  based  on 
experimental  evidence  (In  general,  the  radiative  transfer  theory 

is  applicable  when  the  particle  density  is  approximately  less  than 
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IX  in  volume.). 

III.  Future  Research  Areas 

Having  developed  our  capability  in  solving  the  vector  transport  equation  for 
nonspherical  particles,  we  are  now  in  a  position  to  expand  our  horizon  to  perforin 
research  in  the  following  areas: 

•  Further  development  of  general  numerical  codes  of  vector  radiative 
transfer  solutions. 

0  Beam  propagation  in  an  ensemble  of  spherical  or  nonspherical  particles 
where  multiple  scattering  effects  are  important. 

•  Wave  propagation  in  a  medium  with  non-uniform  density  distribution 
of  particles. 

•  Pulse  propagation  in  a  random  distribution  of  discrete  spherical  or 
nonspherical  scatterers. 

IV.  Listings  of  Papers  and  Presentations 

The  following  is  a  list  of  papers  and  presentations  by  our  group  during  the  current 
ARO  research  program. 

(a)  "Matrix  representations  of  the  vector  radiative  transfer  theory  for  ran¬ 
domly  distributed  nonspherical  particles"  JOSA  A,  1_,  pp.  359-364  (1984). 

(b)  "First  order  multiple  scattering  theory  for  nonspherical  particles",  to 
appear  in  Appl.  Opt.  (1984) 

(c)  "Multiple  scattering  calculations  for  nonspherical  particles"  based  on 
the  vector  radiative  transfer  theory"  to  appear  in  Radio  Science  (1984). 

(d)  "Multiple  scattering  by  nonspherical  particles"  presented  at  the  1982, 

1983  and  1984  CLS  Scientific  Conference  on  Obscuration  and  Aerosol 
Research,  Aberdeen  Proving  Group  (U.  S.  Army),  Maryland. 


(e)  "Vector  radiative  transfer  theory"  presented  at  the  1984  URSI 
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